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Theory of the Single-Material, Helicoidal Fiber 

By J. A. ARNAUD 
(Manuscript received March 29, 1974) 

The theory of propagation in a new single-material, single-mode, optical 
fiber is given. The modes are of the ivhispering-gallery type, with the 
propagation taking place along helicoidal paths close to the boundary of a 
cylindrical dielectric rod. The beams are confined in the azimuthal direc- 
tion in helicoidal ridges. It is shown that single-mode, low-loss operation 
is possible if the helix period is of the order of the rod cross-section area 
divided by the wavelength and the ridge area is of the order of 1 percent of 
the rod cross-section area for two channels. The rod is supported by heli- 
coidal wings that play a role in the mode-selection mechanism. 

I. INTRODUCTION 

The best-known single-mode optical fiber is the clad fiber. If the 
difference in refractive index between core and cladding is small, 
single-mode propagation can be achieved for core diameters that are 
large compared with the wavelength. It is, however, desirable to use 
just one material, such as quartz, that exhibits low impurity and scat- 
tering losses. In a previous work, 1 we indicated that single-mode propa- 
gation could be achieved in a single-material configuration that we 
called a "helicoidal fiber." Figure 1 represents a more recent version 
of this type of fiber. 

To explain the mechanism of operation, let us consider first a cy- 
lindrical dielectric rod with radius a = p. The refractive index of the rod 
is perhaps n = 1.45 (quartz), and the surrounding medium is air. Waves 
are guided along the rod boundary as shown in Fig. 2a. These so-called 
"whispering-gallery modes" can be represented by rays repeatedly 
reflected from the boundary because of total reflection. In the interior 
of the rod, the modes are described by Bessel functions J v (kr) 
X exp (iv<t>), where v is a large integer and kr is a large number of the 
order of v. Because kr and v are both large and comparable to one 
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Fig. 1 — Open view of the single-material, helicoidal fiber for two optical channels. 
The two optical beams propagate in the two ridges shown, with areas t X 2d. The 
helicoidal motion is essential to maintain confinement. High-order modes are not 
confined. They radiate away to the envelope through the wings (only part of one is 
shown). The ratio period /diameter is much larger than that shown in the figure. 

another, the Bessel functions can be approximated by Airy functions. 
The field is oscillatory from the rod boundary down to a slightly smaller 
radius r c called the caustic (or turning point) radius. For radii smaller 
than r c , the field decays exponentially. Thus, the field of whispering- 
gallery modes clings tightly to the rod boundary. The distance between 
the caustic and the boundary, which defines in some sense the "thick- 
ness" of the mode, is for the fundamental mode of the order of (X 2 a)*, 
with X the wavelength in the medium and a the rod radius. For ex- 
ample, if X = 1 fxm. and a is equal to 8 mm, the fundamental mode 
thickness is of the order of 20 /jm. It increases with the mode num- 
ber m, approximating as w J . As m increases, the phase velocity in- 
creases too. 

These whispering-gallery modes can be generalized to take into 
account a motion along the rod axis z. The combined rotation and axial 
motion results in a helicoidal path that can be understood from simple 
ray-optics considerations. The only significant difference from the 
previous case is that the radius a of the rod should be replaced, in the 
expression for the mode thickness, by the helix radius of curvature, 
p = a/sin 2 8, where denotes the angle that the helix makes with the 
rod axis. For example, if a is 80 /mi and 6 = 0.1 radian, the mode 
thickness is the same as in the previous example, where a was assumed 
to be 8000 fj.m. If is equal to zero, there is of course no confinement 
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Fig. 2 — (a) Whispering-gallery modes clinging to a circular boundary with effective 
radius p( = a/sin 2 0). The field (m = 0, 1, • • •) is described exactly by Bessel functions 
and approximately by Airy functions, (b) Cross section in a local r, £ plane. 



near the rod boundary. Observation of helicoidal rays in optics has 
been reported. 2 

Let us now assume that we have selected a convenient value for 0, 
perhaps 6 = 2.5°, and that we wish to define one or more channels in 
the azimuthal direction. This can be achieved with helicoidal sepa- 
rators, 1 or ridges, as shown in Fig. 1, that follow the path of the desired 
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whispering-gallery modes. It is clear, intuitively, that the optical power 
will tend to remain in the ridges. As the mode number increases, either 
in the azimuthal or radial direction, the modes occupy larger and 
larger volumes and eventually "spill out" of the ridges. Because there 
is maximum confinement in the ridge when 6 = ir/2 and no confine- 
ment at all when 6 = 0, it is plausible that only a single mode remains 
confined in the ridge for a proper choice of 0. The higher-order modes 
radiate away from the ridge along the boundary. They can be ab- 
sorbed easily without degrading the fundamental mode. In this paper, 
we justify the above intuitive arguments and show that strong dis- 
crimination against unwanted modes can indeed be obtained. 

The first step of the calculation is to obtain the propagation con- 
stants of whispering-gallery modes in circular cylinders in a con- 
venient form. This is done in Section II. In Section III, we investigate 
the case of helicoidal boundaries in the local mode approximation and 
obtain the design parameters. In Section IV, the case of small ridges 
is investigated with the help of a new perturbation method. 

The single-material helicoidal fiber discussed in this paper can be 
compared to the single-material ridge guide recently demonstrated 3 
and analyzed. 45 These two single-material fibers have features in 
common. The mode-selection mechanism rests on similar general 
principles. It can be ascribed to a coupling between ridges carrying 
trapped modes and two-dimensional substrates carrying radiation 
modes. 6 In the case of the ridge guide, the slab constitutes the two- 
dimensional substrate needed to ensure single-mode propagation. In 
the case of the helicoidal fiber, the dielectric rod itself can be con- 
sidered a two-dimensional mode sink because the whispering-gallery 
modes that it guides have a restricted thickness in the radial direction, 
as we have discussed before. In both cases, a good discrimination 
against high-order modes should in principle be obtained by increasing 
the distance between the absorbing elements and the ridges, because 
these elements are coupled through the radiation field rather than 
through evanescent waves. 

The theory given in this paper is applicable to purely metallic heli- 
coidal waveguides as well as to dielectric waveguides. The metallic 
helicoidal waveguide is attractive as a low-loss, multichannel, single- 
mode system for long-distance microwave communication. It can be 
compared to the groove guide 7 shown in Fig. 4c. The metallic helicoidal 
waveguide has the advantage that TEM modes are absent. In the 
groove guide, any lack of symmetry between the two plates introduces 
a large loss through coupling to the (slower) TEM modes. This is, in 
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fact, also the reason for the superiority of the dielectric ridge guide 4 
over the metallic groove guide. 7 

II. PROPAGATION OF WHISPERING-GALLERY MODES IN CYLINDRICAL 
SURFACES 

Let us consider a circular dielectric cylinder with radius a. For E 
modes, the </> and z components of the electric field have the form 

E(r, </», z) = J t (ur/a) exp [i(v+ + ft«aQ, (1) 

where 

U 2 = ( fc 2 _ jfi a i m ( 2 ) 

Assuming that the discontinuity in refractive index is sufficiently 
large, a condition well satisfied for quartz rods in air, the boundary 
condition at r = a is 

J,(u) = 0, (3) 

because, for the type of mode considered, the field tends to vanish at 
the boundary. The zeros of J v are denoted u m (v), m = 0, 1, 2, • ■ ■. 
We introduce new coordinates £, f in place of y, z (see Fig. 1) 

£ = r cos 6 <f> — sin 6 z, ... 

(4) 
f = cos 9 z + r sin 6 <t>, 

where 

6 m tan" 1 (2*r/p), (5) 

the quantity p, for "period," being for the moment an arbitrary con- 
stant. The wave numbers r ? , r r in the new coordinate system are 
related to v, k z by 

v = r cos 6 T f + ?• sin 6 T r , (6a) 

k z = -sin r £ + cos I>. (6b) 



They are such that 



I\£ + T f f = v<f> + k z z. (6c) 



The characteristic equations, (2) and (3), are now written, using 
eqs. (6), 

(-sin r £ + cos I» 2 a 2 + u 2 n (r cos r* + r sin I» = k 2 a 2 . (7) 

Equation (7) provides us with the desired relation between r f and r { . 
We wish to simplify this relation. Because we are interested in whisper- 
ing-gallery modes corresponding to large values of v, we can use the 

SINGLE-MATERIAL, HELICOIDAL FIBER 1647 



Table I — Values of b parameter in eq. (10) 



m 


b m 


[37r(2m + 3/4)/2»4]I (J.W.K.B.) 



1 
2 
3 


1.85575 
3.24461 
4.38167 
5.387 


1.841 
3.239 
4.379 
5.385 



following approximation for u m (v) 8 



u m (v) = v + b m v\ 



(8) 



where b m is given in Table I. In the second column in Table I, the 
J.W.K.B. approximation for b m obtained from simple ray optics con- 
siderations is given. As we can see, the error does not exceed 1 percent 
even for small m. 

We note further that a is not very different from r. Thus, we set 
a = r -f- t, t <3C r. Because we are considering waves that do not depart 
very much from the reference helicoidal path, T f is very close to k, and 
the transverse wave number r { is small compared with iy Neglecting 
products of small quantities, (7) becomes 



r» a r? + rf = fc 2 [i + 2t/ P - 2b m {k P )-*} i 



(9) 



where we have introduced the reference helix curvature p = a/sin 2 6. 
The term 2t/p expresses the fact that, at the reference radius r, the 
phase velocity is smaller than at the boundary with radius a. The term 
2b m (kp)~ i results from the radial variation of the field. The larger the 
radial mode number m, the smaller the tangential wave number r. 
Note that the system is approximately isotropic. 

III. HELICOIDAL BOUNDARY 

In the previous section we have assumed that the boundary is a 
circular cylinder with radius a. We now assume that a is a function of 
£, but that it remains independent of f . By letting a vary with £, we 
generate a helicoidal surface. Azimuthal confinement of the whispering- 
gallery beams can be expected for various well-shaped profiles a(£). 
For simplicity, we assume here that a(£) = a, a constant, for — d < £ 
< d, and a = r, where r denotes the reference radius, anywhere else 
in the period. A slightly tapered transition region is assumed. Mode 
mixing can therefore be neglected in the evaluation of the propagation 
constants of the modes m = 0, 1, 2, • ■ • . A small amount of mode mix- 
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ing is nevertheless needed for this mode-selection mechanism to 
operate. 4-7 

From (9), the wave numbers for the fundamental mode (m = 0) 
and first-order mode (m = 1) in the ridge (unprimed number) and 
outside the ridge (primed number) are, respectively, 

Tl = /c 2 [l + 21/ p - 26 (fcp)-»], (10a) 

Tl> = /c 2 [l - 2bo(fcp)-*} (10b) 

r? - fc 2 [l + 2t/p - 26i(fcp)-»], (10c) 

r 2 , = /c 2 [l - 26i(fcp)-»]. (lOd) 

The axial wave numbers r ro and r f i for the two modes and 1 are 
now obtained using the standard dielectric slab theory. If we normalize 
the axial wave number r f0 by defining 

K> m (r 2 - rgo/(rg - rfr) (li) 

and introduce the V parameter 

F 2 = (rg - Tl)d\ (12) 

we have, for the modes m = 0, an explicit relation between V and K, 

V = {tan- 1 [X(l - K 2 )-^ + wtt/2}(1 - K*)-\ (13) 

where n = 0, 2, 4, • • • correspond to even modes and n = 1,3, • • • to 
odd modes.* A similar relation holds for the modes m = 1, n = 0, 
1, 2, ■■. 

The axial wave numbers r fmn of these various modes m,n are plotted 
in Fig. 3 as functions of the ridge width 2d, for X = 1 urn, n = 1.45 
(quartz), a rod radius a = 50 nm, and a ridge height t = 3.5 /im. We 
have chosen 6 = 2.5°, corresponding to a helix radius of curvature 
p = 25 mm. Modes whose axial wave number is less than IY (the wave 
number of the fundamental mode outside the ridge) suffer radiation 
losses. 

This figure clearly shows that only one mode (m = 0, n = 0) is free 
of radiation loss if 2d is less than 14 pm. For 2d = 14 pm, the field of 
the fundamental mode decays in azimuth by a factor of 1/e at a 
distance 

£ = (I> 2 o - rgO-i = 6.3 pm (14) 

on either side of the ridge. For two channels, the "wings" holding the 



* The mode number n should not be confused with the refractive index n. 
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Fig. 3 — Axial wave numbers r fmn for various m and n modes (r and £ coordinates). 
Radiation loss is suffered whenever Tf < IV. The radiation zone is shaded. This 
figure shows that single-mode propagation is possible if the ridge width 2d is less 
than 14 /un (for X = 1 ^m, n = 1.45, t = 3.5 ftm, o = 50 nm, and = 2.5°). 

rod in Fig. 1 are located at a distance Ta/2 = 80 fxra from the corruga- 
tion. At that distance, the field has decayed by a factor of more than 
10 5 . The fundamental mode therefore suffers negligible radiation loss 
(bending losses are not considered here). If we set the condition that 
£o be tg of 7ra/4, we obtain the approximate condition ^ \/a for 
single-mode low-loss operation. More detailed relations are given at 
the end of Section IV. 

The local mode theory used in this section is expected to be appli- 
cable when the ridge width 2d is large compared to the ridge height, t, 
and large compared to the wavelength, X. In the next section, we find 
that a perturbation method applicable to small ridge areas (2 td) leads 
to almost identical conclusions. 

IV. LINE PERTURBATION OF SURFACE WAVES 

We give a general theory of the trapping of surface waves by rods 
of small cross section. This theory is then applied to helicoidal ridges 
of small cross section. 

Let us consider an isotropic surface, perhaps an inductive surface, 
supporting a plane wave with wave number k. Let us introduce a di- 
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Fig. 4 — Various methods of introducing a line perturbation on a surface wave, 
(a) Reactive (e.g., corrugated) surface perturbed by a dielectric rod. (b) TEM waves 
perturbed by a dielectric slab, (c) "Groove guide." 7 (d) Ridge guide considered in 
this paper. The torsion of the helicoidal motion is not essential. The radius of curva- 
ture p of the helix is the important parameter that determines mode selection. 

electric rod of very small cross section parallel to this surface, some 
distance away from it (Fig. 4a). Because the power carried by the plane 
wave is infinite, a straightforward application of the conventional 
perturbation method does not give any meaningful result. Therefore, 
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we shall proceed the other way around. We start from the perturbed 
state and assume that we know the propagation constant h > k. The 
wave is in that case confined transversely with a decay rate (h 2 — A; 2 ) 1 . 
We now "peel off" the rod and evaluate the successive perturbations 
until the perturbation resulting from the rod vanishes. By specifying 
that h — > k in that limit, the initial value of h is obtained. 

For simplicity, this method is first explained for the case where 
e £^ € , the perturbed field being of the order of the unperturbed field. 
Let a be a parameter such that a = corresponds to the absence of the 
rod and a = 1 corresponds to the presence of the rod. Furthermore, a 
is so chosen that, in the perturbation formula 

dh/dff = a(h 2 - k 2 )i, (15) 

a is a constant. This can be done because we have factored out a term 
(h? — k 2 ) i inversely proportional to the power carried by the mode. 
(The distortion of the field in the close neighborhood of the rod does 
not contribute significantly to the total power, because of the large 
transverse extent of the field.) The ratio <r is essentially the ratio of the 
present cross-section area of the perturbating rod to its original cross- 
section area. Integrating eq. (15) from a = 1 to a = 0, we obtain 

f ' (h 2 - k 2 )~Uh = a (16a) 

or 

h = k cos a PS fc(l - a 2 /2). (16b) 

To clarify the significance of this result, let it be applied to a con- 
figuration where the exact solution is known. Consider two parallel 
perfectly conductive plates with spacing D carrying TEM modes, as 
shown in Fig. 4b. If we introduce a dielectric slab with e PS « and 
width 2d, we obtain the so-called "H-guide" configuration proposed by 
Tischer. 9 (Note, however, that we consider here the H modes rather 
than the low-loss modes.) The parameter a = y/d, where y is shown 
in Fig. 4b, clearly satisfies the requirements set up above. The con- 
ventional perturbation formula (see, for example, Ref. 5, Part II, eq. 
(21), with E p W E, H p « H, E* = E, W = -H) is 

Ah = $a /"(« - e )E 2 dS / JE X H-dS. (17) 

For our case, we obtain, taking into account the exp [— (h 2 — k 2 )* | y | ] 
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dependence of the field on y, 

dh/da = (no/e )Hh 2 - fc 2 )*(« - «oM = a(h? - k 2 )K (18) 

Thus the constant a is, from (18), 

a = (n 2 - l)fcd, (19) 

where e/e = ro 2 . Application of (16b) now gives the perturbed wave 
number h 

h = fc[l + |(n> - l) 2 fc 2 d 2 ]. (20) 

The exact solution to the problem, for small (n 2 — l)fcd, is well known 
[see, for example, Ref. 5, Part II, footnote after eq. (11)]. We have, 
with the approximation tan [(?i 2 — l)'fcd)] « (n 2 — l)*fcd, 

/,2 _ ^2 = („2 _ 1)2^4^2 (21) 

Equation (21) coincides with our perturbation result, (20), because 
h ^5 k. Having satisfied ourselves with the validity of our perturbation 
technique, we apply it to a small wall perturbation. We assume that 
the case of quartz in air is the same as the case of a metallic boundary, 
except for the wavelength \/n replacing X. 

For a wall perturbation with cross-section area s, the perturbation 
formula is 

Mi = iunoH 2 s I /*E X H-dS, (22) 

if the electric field is equal to zero. Note that h increases if the volume is 
increased, e.g., if we introduce corrugations in the wall. 

For H waves uniform along the ?/-axis (see Figs. 4c and 4d) (E u = E), 
we have, from Maxwell's equations, 

H x = - (k/ mo )E, (23a) 

H 2 = (itait )- l dE/dx. (23b) 

Substituting in (22) we obtain the perturbation 

Ah = (s/2k)(dE/dx) 2 (h 2 - fc 2 )* / f&dx. (24) 

Defining a as x/t (see Figs. 4c or 4d), the constant a defined in (15) is 
found to be 

a = (s/2k)(dE/dxy/ JEHx, (25) 
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the derivative being evaluated at the first zero of E(x) for the mode 
m = 0, the second for the mode m = 1, and so on. 
For whispering-gallery modes, we have 

E = Ai(0, (26a) 

where Ai( ) denotes the Airy function and 

t m (2k*/p)*x - (2k 2 /p)-*(k 2 - h 2 ), (26b) 

p being the boundary radius. Substituting in eq. (25), we obtain 

a = 2f m ktd/ P) (27) 

where f m is a numerical factor 

f m = (dAi/dt)f =lm / f tm Ai 2 (t)dt, (28) 



^o, h, ■ ■ ■ ,t m , ■ ■ ■ being the zeros of Ai(t). By numerical integration, we 
find 

/„ = 0.981 • • • (29a) 

/i = 0.955 •••. (29b) 

Thus, the change in propagation constant resulting from a wall 
deformation of area s = 2td is, from (27), 

h m -k = \f 2 m ¥sy P \ (30) 

As long as the perturbation is small, h\ remains smaller than IV and 
only the mode m = is free of radiation loss. For a sufficiently large 
perturbation, however, hi may exceed IV. Then the modes m = and 
m = 1 are both free of radiation loss, that is, the system is no longer 
single-mode. The condition for the system to be single-mode is therefore 

h m - T v < IV - IV, 

or 

/£*V/p* < 2k(b - 6 1 )(/cp)- J - (31) 

Tc, IV, and the constants b , &i were defined in Section III. The above 
condition can be written, using the values in Table I for b 0> b h 

fc 2 s < 1.74(fcp)«. (32) 

If we take the limit d — > in the expressions given in Section III, we 
obtain instead 

k 2 s < 1.68(*p)», (33a) 
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which is very close to the perturbation result, (32). Thus, the local mode 
approach and the perturbation approach agree closely, not only in 
form, but also numerically. Because p £zl a/6 2 , 9 being the angle that 
the helix makes with the z axis, and 8 ^ 2-wa/p, p being the helix 
period, condition (33) for single-mode operation can be rewritten 

s < 0.012 (xy/o)». (33b) 

The extent £ in azimuth of the fundamental mode, defined as the 
1/e point of the field, £ = (hi - fc 2 )"* PS (2A0 _i (^o - fc)~*, is, from 
(30) with m = 0, given by 

fc£o « p/ks. (34a) 

If we specify that the fundamental mode field has decayed by a factor 
of 10 5 at the "wings," located, for two channels (see Fig. 1), a distance 
to/2 away from the ridge, we must have £ = (ira/2)/ll.5. Introducing 
the helix period p, this condition for the fundamental mode to have 
small radiation loss can be written 

s > 0.005 (Xp/a) 2 . (34b) 

The condition for a single mode to propagate, (33b), and for the fun- 
damental mode to have small radiation losses, (34b), are consistent if 

p < 5a 2 /X. (35) 

For example, if a = 50 /xm, X = (1/1.45) /*m, according to eq. (35), 
the helix period, p, must be smaller than 18 mm. A period of 10 mm, 
for instance, would be quite adequate. Note that for such rather long 
periods the optical path is not significantly increased by the circular 
motion. Radiation into free space is negligible, as long as the medium 
surrounding the ridge is air. For mechanical reasons, however, we may 
want to use a material with lower index. In that case, radiation into 
the surrounding medium may be a limitation. 

V. CONCLUSION 

The single-material helicoidal fiber proposed earlier by the author 
has been shown to support only one mode, with low radiation loss, 
provided the following two conditions are satisfied : 

Helix period tt rod cross-section area/ wavelength. 
Ridge area £zl rod cross-section area/70. 

More detailed calculations, similar in spirit to the ones given in 
Ref. 5, would be necessary to specify the magnitude of the radiation 
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losses and the exact value of the mode discrimination. The helicoidal 
fiber, like any single-mode fiber with large mode cross section, may be 
sensitive to bending losses. The bending loss is therefore another key 
point that needs to be investigated. 
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